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The theory of the solar/stellar activity cycles is presented, based on the mean-field concept in magnetohydrodynamics. A new approach 
to the formulation of the electromotive force as well as the theory of differential rotation and meridional circulation is described for 
use in dynamo theory. Activity cycles of dynamos in the overshoot layer (BL-dynamo) and distributed dynamos are compared, with 
the latter including the influence of meridional flow. The overshoot layer dynamo is able to reproduce the solar cycle periods and the 
butterfly diagram only if a = in the convection zone. The problems of too many magnetic belts and too short cycle times emerge if 
the overshoot layer is too thin. The distributed dynamo including meridional flows with a magnetic Reynolds number Rm 20 (low 
magnetic Prandtl number) reproduces the observed butterfly diagram even with a positive dynamo-a in the bulk of the convection 
zone. 

The nonlinear feedback of strong magnetic fields on differential rotation in the mean-field conservation law of angular momentum leads 
to grand minima in the cyclic activity similar to those observed. The 2D model described here contains the large-scale interactions 
as well as the small-scale feedback of magnetic fields on differential rotation and induction in terms of a mean-field formulation (A- 
quenching, o-quenching). Grand minima may also occur if a dynamo occasionally falls below its critical eigenvalue. We expressed 
this idea by an on-off a function which is non-zero only in a certain range of magnetic fields near the equipartition value. We never 
found any indication that the dynamo collapses by this effect after it had once been excited. 

The full quenching of turbulence by strong magnetic fields in terms of reduced induction (a) and reduced turbulent diffusivity (tjt) is 
studied with a ID model. The full quenching results in a stronger dependence of cycle period on dynamo number compared with the 
model with o-quenching alone giving a very weak cycle period dependence. 

Also the temporal fluctuations of a and from a random-vortex simulation were applied to a dynamo model. Then the low 'quality' 
of the solar cycle can be explained with a relatively small number of giant cells acting as the dynamo-active turbulence. The simulation 
contains the transition from almost regular magnetic oscillations (many voilices) to a more or less chaotic time series (very few 
vortices). 



1 Introduction 

Explaining the characteristic period of the quasicyclic activity oscillations of stars with the Sun 
included is one of the challenges for stellar physics. The main period is an essential property of 
the dynamo mechanism. Solar dynamo theory is reviewed here in the special context of the cycle- 
time problem. The parameters of the convection zone turbulence do not easily provide us with the 
22-year time scale for the solar dynamo. Even for the boundary layer dynamos, it is only possible 
if a 'dilution factor' in the turbulent electromotive force smaller than unity is introduced which 
parameterizes the intermittent character of the MHD-turbulence. 

A notable number of interesting phenomena have been investigated in the search for the solution of 
this problem: flux-tube dynamics, magnetic quenching, parity breaking, and chaos. Nevertheless, 
even the simplest observation - the solar cycle period of 22 years - is hard to explain (cf. DeLuca 
and Oilman, 1991; Stix, 1991; Oilman, 1992; Levy, 1992; Schmitt, 1993; Brandenburg, 1994a; 
Weiss, 1994). How can we understand the existence of the large ratio of the mean cycle period 
and the correlation time of the turbulence? Three main observations are basic in this respect: 

• There is a factor of about 300 between the solar cycle time and the Sun's rotation period. 

• This finding is confirmed by stellar observations (Fig. [T]). 

• The convective turnover time near the base of the convection zone is very similar to the solar 
rotation period. 
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Figure 1: Stellar cycles: the ratio of the cycle time to the rotation period after Saar and Branden- 
burg (1999). The solid line fits the observations for all stars and leads to the weak proportionality 
Wcyc oc i^^^'^ while a previous fit (dashed line) provided Wcyc oc Q^^''^ (Brandenburg et ai, 1998). 

The problem of the large observed ratio of cycle and correlation times, 

^ > 102, (1) 
''"corr 

constitutes the primary concern of dynamo models. In a thick convection shell this number reflects 
(the square of) the ratio of the stellar radius to the correlation length and numbers of the order 100 
in ([I]) are possible. For the thin boundary layer dynamo, however, the problem becomes more 
dramatic and is in need of an extra hypothesis. 

The activity period of the Sun varies strikingly about its average from one cycle to another. Only a 
nonlinear theory will be able to explain the non-sinusoidal (chaotic or not) character of the activity 
cycle (Fig. ^). A linear theory is only concerned with the mean value of the oscillation frequency. 



2 Basic Theory 

Kinematic dynamo theory utilizes only one equation to advance the mean magnetic field in time, 
i.e. 

^ = rot((u)x(B)+5). (2) 

Here only a non-uniform rotation will be imposed on the mean flow, (u) ; any meridional flow shall 
be introduced later. The turbulent electromotive force (EMF), £ = (u' x B'), contains induction 
Uij and dissipation r/jj^, i.e. 

£i = aij{Bj) + r]ijk{Bj)^k + ... ■ (3) 

Both tensors are pseudo-tensors. While for rjijk an elementary isotropic pseudo-tensor exists 
C^ijk"), the same is not true for . An odd number of il's is, therefore, required for the a-tensor. 



2 





T (YEARS) 



Figure 2: LEFT: The distribution of the solar cycle length does not approach a Dirac function, 
the 'quality' of the cycle only gives values of about 5. RIGHT: The wavelet spectrum of the 
sunspot-number time series shows two peaks for both 10 yr and 100 yr (Frick et ai, 1997a). 



which is only possible with an odd number of another preferred direction, (say) g. The a-effect 
can thus only exist in stratified, rotating turbulence. The first formula reflecting this situation. 



a 



•"corr ^' 
Ca— COS I 



(4) 



was given by Krause (1967) with being the angular velocity of the basic rotation, 6 the colat- 
itude, and Hp the density scale height. Evidently, a is a complicated effect, where the effective 
a might really be very small; the unknown factor Cq in (^ may be much smaller than unity. 
The strength of this effect was computed in recent analytical and numerical simulations for both 
convectively unstable as well as stable stratifications. While Brandenburg et al. (1990) worked 
with a box heated from below, Brandenburg and Schmitt (1998) considered magnetic buoyancy 
in the transition region between the radiative solar core and the convection zone, Brandenburg 
(2000) probed the Balbus-Hawley instability for dynamo-a production. Ferriere (1993), Kaisig 
et al. (1993), and Ziegler et al. (1996) used random supernova explosions to drive the galactic 
turbulence. The magnitudes of the a-effects do not reach the given estimate in these cases: the 
dimensionless factor Cq, seems really to be much smaller than unity. 



2.1 Simple Dynamos 
2.1.1 Dynamo waves 

An illuminating example of kinematic dynamo theory for the cycle period is the dynamo wave 
solution. In plane Cartesian geometry there is a mean magnetic field subject to a (strong) shear 
flow and an a-effect - all quantities vary only in the z-direction with a given wave number (Parker, 
1975). The amplitude equations can then be written as 

A + A = Cc,B, B + B = iCnA, (5) 
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with A representing the poloidal magnetic field component and B the toroidal component. Ca is 
the normalized a and Cq is the normalized shear. The eigenfrequency of the equation system is 
the complex number 




A marginal solution of the magnetic field is found for a 'dynamo number' V = CaCn = 2. In 
that case the field is not steady but oscillates with the (dimensionless) frequency cjcyc = 1- 
The (re-normalized) cycle period is thus given simply by combining the eddy diffusivity and the 
wave number, i.e. LJcyc — ryx^:^ (which is simply the skin-effect relation). With a mixing-length 
expression for the eddy diffusivity, r/x = Cr,/corr/''"corr, one finds the basic relation between the 
cycle time and the correlation time as 

''"eye 1 

''"corr ^iTTCfj 

Note that < 0.3. The ratio between the global scale, R, and the correlation length of the 
turbulence, /corr. determines the cycle time. As this ratio has a minimum value of 10, it is thus no 
problem to reach a factor of 100 between the cycle period and the correlation timcQ The rotation 
period does not enter - it only influences the cycle period in the nonlinear regime (for dynamo 
numbers exceeding 2). For such numbers the frequency increases - and the cycle time becomes 
shorter. 



R 

^eorr 



(7) 



2.1.2 Shell dynamo 

The simplest assumptions about the a-effect and the differential rotation are used in the shell 
dynamo, i.e. 

a = oq cos 9, $7 = const. + Jl'x (8) 

for Xj < X < 1 (Roberts, 1972; Roberts and Stix, 1972). Positive W means radial super-rotation 
and negative Q' means sub-rotation. The solution with the lowest eigenvalue for the sub-rotation 
case in a thick shell is a solution with dipolar symmetry, which is, in the vicinity of the bifurcation 
point, the only stable one (Krause and Meinel, 1988). 

The cycle time grows with the linear thickness D of the shell. Compared with the dynamo wave, 
we find a reduction of (^) by the normalized thickness, d < 0.5. As long as convection zones of 
main-sequence stars are considered, that is not too dramatic, but what about rather thin layers like 
the solar overshoot region? As the cycle period is found to vary linearly with D in similar shell 
models, 

RD 

r ~ 0.26 , (9) 

the cycle time must become dramatically short in very thin boundary layers. 

'On the other hand, in order to reproduce a factor of order 10^ the correlation length must not be too small. 
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2.1.3 The dynamo dilemma 

The sign of a has been presumed as positive in the northern hemisphere and the differential rota- 
tion has been considered unknown (Steenbeck and Krause, 1969) from the early years of dynamo 
theory. Only sub-rotation - with angular velocity increasing inwards - could then produce a suit- 
able butterfly diagram from the latitudinal migration of the toroidal field. In the other case, e.g. 
where the surface rotation law is applied in the entire convection zone, one cannot reproduce the 
migration of the toroidal field towards the equator (Kohler, 1973). 

The correct field migration according to the sunspot butterfly diagram is generally produced by 
a negative dynamo number only, which may be either due to positive a in the northern hemi- 
sphere and sub-rotation, or negative a and super-rotation. For all dynamos with super-rotation, 
however, the phase relation between the radial field component and the toroidal field component 
is {Br){B^) > 0, again disagreeing with the observations (Stix, 1976). The prediction of dynamo 
theory for the solar differential rotation was thus a clear sub-rotation, dQ/dr < 0. This prediction 
contrasts with the results of the helioseismology which revealed sub-rotation only near the poles, 
whereas near the solar equator there is a super-rotation at the base of the convection zone (Fig. ||). 
Agreement with the butterfly diagram shall thus be satisfied by a negative a which can be only 
explained for the overshoot region. But the phase relation of the field components does then 
not agree with the known one (Yoshimura, 1976; Parker, 1987; Schlichenmaier and Stix, 1995). 
With the incorrect phase relation it is even problematic to derive the observed properties of the 
torsional oscillations (Howard and LaBonte, 1980), which are suggested to arise as a result of the 
backreaction of the mean-field Lorentz force (Schiissler, 1981; Yoshimura, 1981; Riidiger et at, 
1986; Kxiker et at, 1996). 

In the light of recent developments for dynamos with meridional flow included, there could easily 
exist a solution of this dynamo dilemma in quite an unexpected way (see Section 4). 

2.2 Differential Rotation 

The theory of the 'maintenance' of differential rotation in convective stellar envelopes might be an 
instructive detour from dynamo theory. Differential rotation is also turbulence-induced but without 
the complications due to the magnetic fields. It is certainly unrealistic to expect a solution of the 
complicated problem of the solar dynamo without an understanding of mean-field hydrodynamics. 
There is even no hope for the stellar dynamo concept if the internal stellar rotation law cannot be 
predicted. 

The main observational features of the solar differential rotation are 

• surface equatorial acceleration of about 30%, 

• strong polar sub-rotation and weak equatorial super-rotation, 

• reduced equator-pole difference in Q at the lower convection-zone boundary. 

The characteristic Taylor-Proudman structure in the equatorial region and the characteristic disk- 
like structure in the polar region are comprised by the results. In the search for stellar surface 
differential rotation, chromospheric activity has been monitored for more than a decade. Sur- 
prisingly enough, there is not yet a very clear picture. For example, the rotation pattern of the 
solar- type star HD 114710 might easily be reversed compared with that of the Sun (Donahue and 
Baliunas, 1992).0 

^if the same butterfly diagram is applied. 
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Figure 3: LEFT: The internal rotation of the Sun after the inversion of SOHO data (Kosovichev 
et al, 1997). RIGHT: Theory of the solar internal rotation by Kitchatinov and Riidiger (1995). 
The rotation frequency is given for the equator (solid), mid-latitudes (dashed) and poles (dashed- 
dotted). 



In close correspondence to dynamo theory we develop the theory of differential rotation in a mean- 
field formulation starting from the conservation of angular momentum, 

— (pr^ sin2 e^l) + — {pr sin OQi^) = 0, (10) 

where the Reynolds stress Qj,^ derived from the correlation tensor 

= K(x,i)^(x,0) (11) 
corresponds to the EMF in mean-field electrodynamics. 

The correlation tensor involves both dissipation ('eddy viscosity') as well as 'induction' (A-effect): 

Qij = J^ijk^k — -^ijkl^k,!- (12) 

Both effects are represented by tensors and must be computed carefully. For anisotropic and 
rotating turbulence the zonal fluxes of angular momentum can be written as 

dn , , / (90 dn\ 

Qrd) = —i'±rsmd— (z/|| — z/x)sin0cos0 rcos^-- sin^-— M- 

or " \ or oO J 

+ +sin2^yW) Osin0, (13) 

Qes = -'^±sin6— - - i'±) sin^ ( sin 6'— - r cos0— ) -|- i^r/f O sin^ 6' cos 6* (14) 

at' " V od or J 

(see Kitchatinov, 1986; Dumey, 1989; Riidiger, 1989). All coefficients are found to be strongly 
dependent on the CorioUs number 

n* = 2rcorr ^ (15) 

with Tcorr as the convective turnover time (Oilman, 1992; Kitchatinov and Riidiger, 1993). More- 
over, the most important terms of the A-effect (y^'^\ H^^^) correspond to higher orders of the 
Coriolis number. The Coriolis number exceeds unity almost everywhere in the convection zone 
except the surface layers. That is true for all stars - in this sense all main-sequence stars are rapid 
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Figure 4: LEFT: The turbulence-originated coefficients in the non-diffusive zonal fluxes of angular 
momentum in the solar convection zone. 

(sohd), (dashed), V^^^ + F^^) (dot-dashed). 
RIGHT: Eddy viscosity in units of 10^^ cm^/s (solid) and Coriolis number Q* in the solar convec- 
tion zone. 



rotators. Theories linear in J7 are not appropriate for stellar activity physics. The a-effect will 
thus never run with cos 6 in its first power. 

The Coriolis number is smaller than unity at the top of the convection zone and larger than 
unity at its bottom. At that depth we find minimal eddy viscosities ('i7-quenching') and maximal 
y{i) ^ ^(1) (pjg |) sjj^£,g the latter are known as responsible for pole-equator differences in Vl, 
we can state that the differential rotation is produced in the deeper layers of the convection zone 
where the rotation must be considered as rapid {Q.* < 10). 

The solution of (|l^) with the turbulence quantities in Fig. ^ is given in Kiiker et al. (1993) and 
Kitchatinov and Riidiger (1995) using a mixing-length model by Stix and Skaley (1990). With 
a mixing-length ratio (omlt = 5/3) we find the correct equatorial acceleration of about 30%. 
There is a clear radial sub-rotation (9J7 /9r < 0) below the poles while in mid- latitudes and below 
the equator the rotation is basically rigid (right panel in Fig. ||). In this way the bottom value of 
the pole-equator difference is reduced and the resulting profiles of the internal angular velocity are 
close to the observed ones. Fig. |5| presents the results of an extension of the theory to a sample of 
main-sequence stellar models given in Kitchatinov and Riidiger (1999). Simple scalings like 

^ oc 17^^ (16) 

may be introduced for both the latitudinal as well as the radial differences of the angular velocity. 
The K-exponents prove to be negative and of order —1 (right panel of Fig. So we find that for 
one and the same spectral type the approximation 

50. ~ const. (17) 

should be not too rough. Recent observations of AB Dor (Donati and Cameron, 1997) and PZ Tel 
(Barnes et al., 1999) seem to confirm this surprising and unexpected result where in all 3 cases the 
constant value approaches 0.06 day~^. 



2.3 The EMF for Turbulence 

We now step forward to the turbulent EMF using the same turbulence model as in Section 2.2. We 
are working with the same quasilinear approximation and the same convection zone model. The 
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Figure 5: The rotation differences 6i^/Q for both a solar model and a K5 main-sequence star vs. 
basic rotation after Kitchatinov and Rudiger (1999). LEFT: Normalized equator-pole difference 
of surface rotation, CENTER: The normalized radial differential rotation, RIGHT: The scaling 
exponent k in (|^) for the solar model. In both cases k ~ — 1 is approximated. 



additional equation for the magnetic fluctuations is 

— - 77tAB' = rot(u' X (B)), (18) 

where rjt ~ vt is the small-scale diffusivity. In this configuration we deal with an a-tensor being 
highly anisotropic, even for the most simple case of slow rotation ($7* ^ 1), 

aij = jEijkGk - ai{GQ)5ij - a2iGiQj + GjQi) + ^EijkUk - aiiUft)5ij - a2iUiQj + UjQi) 

(19) 

(Moffatt, 1978; Krause and Radler, 1980; Walder et al, 1980) with the stratification vectors 
G = V log p and U = V log with the rms velocity = \/ (u'^). The 7-terms describe 
advection effects, such as turbulent diamagnetism (Radler, 1968; Kitchatinov and Rudiger, 1992) 
and buoyancy (Kitchatinov and Pipin, 1993) and cover the off-diagonal elements of the a-tensor. 
The diagonal elements are essential for the induction. The mixing-length approximation used in 
Rudiger and Kitchatinov (1993) leads to 

Urr = a + 17, = aee = -a + fi. (20) 

While the most important component q^^ becomes positive (if density stratification dominates), 
the component a^r becomes negative in the northern hemisphere. In contrast to the standard 



formulation, the a2-components in ( |19| ) are dominant as confirmed by numerical simulations by 
Brandenburg et al. (1990). Our quasilinear theory of the a-effect provides a = 8/15 t'^^„u\. 
The general case of an arbitrary rotation rate is not easy to present. The 'overshoot' region, 
however, is of particular interest and is characterized by |U| > |G| and 17* » 1. In cylindrical 
coordinates (s, cf), z) and for very fast rotation one gets 

fcos0 Ifrcos^ v.. 



a = Cn I — ^fr COS 9 — COS b 

-If. sin I 



dr 



(21) 



(Rudiger and Kitchatinov, 1993). Note that 
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Simulation of the a-effect after Brandenburg (1994b) versus latitude and depth. Note the 
values at the bottom of the convection zone (RIGHT) and the maximum for high latitudes 



• all components with index z disappear, 

• the remaining diagonal terms (the a-effect) do not vanish for very rapid rotation, 

• the a-terms are negative (in the northern hemisphere) for outward increase of the turbulence 
intensity (like in the overshoot region), 

• the advection terms tend to vanish for rapid rotation; their relation to the a-effect is given 
by the Coriolis number Vi* . 

It might seem that the a-effect is now in our hands. There are, however, shortcomings beyond the 
use of the quasilinear approximation. The main point is the absence of terms higher than first order 
in G or U, respectively. In particular, all terms with (Gfi)'^ ~ cos'^ Q are missing so that we have 
no information about the a-effect at the poles where for slow rotation a^^ is largest. However, 
this need not to be true for 'rapid' rotation. Simulations of a^^ with its latitudinal profile by 
Brandenburg (1994b) indeed showed that the maximum occurs at mid-latitudes (Fig. |6|). Schmitt 
(1985, 1987) derived a similar profile for a different type of instability. 

A new discussion of the a-effect in shear flows recently emerged (cf. Brandenburg and Schmitt, 
1998). A quasilinear computation of the influence of the differential rotation on the a-effect leads 
to 

2 517 d log (/jut) . . 

« - -^corr^ Sm^ (22) 

for a sphere or 

(Kl d log p 

ds dz 

for a disk (Riidiger and Pipin, 2000). The latter relation for accretion disks with dQ./ds < yields 
negative a-values in the northern hemisphere and positive values in the southern hemisphere (see 



~ _/2 u-g L u. p 



Brandenburg and Donner, 1997). According to (|2lD or (g2|), negative values for the northern a^,^ 
in the solar overshoot region are only achievable with a very steep decrease of the turbulence 
intensity with depth. 

There are very interesting observational issues concerning the current helicity 

Wcurr = if ■ B') (24) 
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with the electric current j = rot B/hq. The current helicity has the same kind of equatorial (anti-) 
symmetry as the dynamo-a. For homogeneous global magnetic fields, the dynamo-a is related to 
the turbulent EMF according to UijBiBj = S ■ B. After Radler and Seehafer (1990) we read this 
equation as a,^^ = S ■ B/B"^, where a^^ is the dominant component of the a-tensor. We are, in 
particular, interested in checking their and Keinigs' (1983) antiphase relation, 

B^ 

= -r? < 0, (25) 



curr 

between a-effect and current helicity. In addition, we have recently shown that indeed the nega 



tivity of the left-hand-side of (25) is preserved for magnetic-driven turbulence fields which do not 



fulfill the conditions for which the Keinigs relation (g5|) originally has been derived (Riidiger et 
al, 2000). 

An increasing number of papers presents observations of the current helicity at the solar surface, 
all showing that it is negative in the northern hemisphere and positive in the southern hemisphere 
(Hale, 1927; Seehafer, 1990; Pevtsov et al, 1995; Abramenko et al, 1996; Bao and Zhang, 1998; 
see Low (1996) for a review). There is thus a strong empirical evidence that the a-effect is positive 
in the bulk of the solar convection zone in the northern hemisphere. 

For a given field of magnetic fluctuations the dynamo-a as well as the kinetic and current he- 
licities, have been computed by Riidiger et al (2000) assuming that the turbulence is subject 
to magnetic buoyancy and global rotation. In particular, the role of magnetic buoyancy appears 
quite important for the generation of a-effect. So far, only the role of density stratification has 
been discussed for both a-effect and A-effect, and their relation to kinetic helicity and turbulence 
anisotropy. If density fluctuations shall be included then it makes no sense here to adopt the 
anelastic approximation div pu' = but one has to work with the mass conservation law in the 
form 

do' 

^ + pdivw' = 0. (26) 
at 

For the turbulent energy equation one can simply adopt the polytropic relation p' = c^cP', where 
Cac is the isothermal speed of sound. The turbulence is thus assumed to be driven by the Lorentz 
force and it is subject to a global rotation. The result is an a-effect of the form 

5 ci, Hop 

Here g is the gravitional acceleration, and the energy of the magnetic fluctuations forcing the 
turbulence is also given. For rigid rotation the a-effect proves thus to be positive in the northern 
hemisphere and negative in the southern hemisphere. Again we do not find any possibility to 
present a negative a-effect in the northern hemisphere. 
The kinetic helicity 

Tikin = • rotu') (28) 

has just the same latitudinal distribution as the a-effect, but the magnetic model does not pro- 
vide the minus sign between the dynamo-a and the kinetic helicity which is characteristic within 
the conventional framework, a oc — Tcorr ^kin- If a rising eddy can expand in density-reduced 
surroundings, then a negative value of the kinetic helicity is expected. The magnetic -buoyancy 
model, however, leads to another result. 
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If the real convection zone turbulence is formed by a mixture of both dynamic-driven and magnetic 
driven turbulence, the kinetic helicity should also be a mixture of positive parts and negative parts 
so that it should be a rather small quantity. The opposite is true for the current helicity which also 
changes its sign at the equator. Again it turns out to be negative in the northern hemisphere and 
positive in the southern hemisphere. The current helicity (due to fluctuations) and the a-effect are 
thus always out of phase. The current helicity at the solar surface might thus be considered a much 



more robust observational feature than the kinetic helicity (28). 

Strong magnetic fields are known to suppress turbulence and the a-effect is also expected to de- 
crease. The field strength at which the a-effect reduces significantly is assumed, in energy, to be 
comparable with the kinetic energy of the velocity fluctuations, and is called the equipartition field 
strength Bcq = (fiopu^)^^ . The magnetic feedback on the turbulence was studied by Riidiger 
& Kitchatinov (1993) using a (5-function spectral distribution of velocity fluctuations. The a^^j,- 
component of the a-tensor acting in an aO-dynamo is found to be suppressed by the magnetic 
field with the a-quenching function 

15 / 4/?2 1 - /32 \ 

with 

13=\{B)\/B,q. (30) 

The quenching function decreases as for strong magnetic fields. The function is similar to 
the heuristic approach (1 + which is artificial but nevertheless frequently used in a-effect 

dynamos. For small fields, ^' runs like 1 — ^/3^, the mentioned heuristic approximation runs like 

1-/32. 

We can also compute the eddy diffusivity using the same turbulence theory as in the previous 
Section. The simplest expression is 

VT = C^TcorrUT (31) 

with < 0.3. This expression, however, fails to include rotation and stratification. Even a non- 
uniform but weak magnetic field induces a turbulent EMF, without rotation and stratification in 
accordance to rjij^ = r]T^ijk- A much stronger 'magnetic quenching' of such an eddy diffusivity 
has been discussed by Vainshtein and Cattaneo (1992) and Brandenburg (1994a) and is still a 
matter of debate. A more conventional theory of the r/x-quenching and its consequences for the 
stellar activity cycle theory is given in Section 6 below. 



3 A Boundary-Layer Dynamo for the Sun 

The spatial location of the dynamo action is unknown unless helioseismology will reveal the exact 
position of the magnetic toroidal belts beneath the solar surface (cf. Dziembowski and Goode, 
1991). There are, however, one or two arguments in favour of locating it deep within or below the 
convection zone: 

• Hale's law of sunspot parities can only be fulfilled if the toroidal magnetic field belts are very 
strong (10^ Gauss, see Moreno-Insertis, 1983; Choudhuri, 1989; Fan et at, 1993; Caligari 
etal, 1995). 
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Figure 7: The depth-dependence of the a-effect with Cq, = 1. Note the different contributions by 
density gradient and turbulence intensity distribution. The a-effect becomes negative only at the 
bottom of the convection zone (Krivodubskij and Schultz, 1993). 



• The dynamo field strength is approximately the equipartition value Bcq = (fiopu^) . 
Using mixing-length theory arguments pu^ ~ const., hence Bcq is increasing inwards with 

(one order of magnitude, < 10^ Gauss). 

• The radial gradient of i7 is maximal below the convection zone. 

As a consequence of these arguments high field amplitudes might be generated only in the layer 
between the convection zone and the radiative interior (Boundary layer (BL), 'tachocline', see van 
Ballegooijen, 1982). On the other hand, if there is some form of 'turbulence' in this layer, it will be 
hard to understand the present-day finite value of lithium in the solar convection zone. The lithium 
burning only starts 40 000 km below the bottom of the convection zone. Any turbulence in this 
domain would lead to a rapid and complete depletion of the lithium in the convection zone, which 
is not observed. Moreover, as shown by Riidiger and Kitchatinov (1997), it is also not easy to 
generate very strong toroidal magnetic fields in the solar tachocline as the Maxwell stress tends to 
deform the rotation profile to more and more smooth functions. In our computation with molecular 
diffusion coefficient values the toroidal field amplitude never did exceed 1 kGauss, independent 
of the radial magnetic field applied. 

Krivodubskij and Schultz (1993) (with the inclusion of the depth-dependence of the Coriolis num- 
ber and using a mixing-length model) derived a profile of the a-effect (Fig. ^ with a magnitude 
of 100 m/s, positive in the convection zone and negative in the overshoot layer. The negativity of a 
there can only be relevant for the dynamo if the bulk of the convection zone is free from a. It has 
been argued that the short rise-times of the flux tubes in the convection zone prevent the formation 
of the a-effect (Spiegel and Weiss, 1980; Schiissler, 1987; Stix, I99I). However, the rise-times 
may be much longer in the overshoot region (Ferriz-Mas and Schiissler, 1993, 1995; van Balle- 
gooijen, 1998). Anyway, as Cq, = 1 is used in Fig. |7|, the given values can only be considered as 
maximal values. 

There is a serious shortcoming of the BL concept. The characteristic scales of the magnetic fields 
are no longer much larger than the scales of the turbulence. The validity of the local formulations 
of the mean-field electrodynamics is not ensured for such thin layers. Nevertheless, there is an 
increasing number of corresponding quantitative models (Choudhuri, 1990; Belvedere et al. , 1991 ; 
Prautzsch, 1993; Markiel and Thomas, 1999). Even in this case a couple of new questions will 
arise. For a demonstration of this puzzling situation a model has been estabUshed by Rudiger and 
Brandenburg (1995) under the following assumptions: 
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• Intermittency (or the weakness of the turbulence in the BL) is introduced by a dilution factor 
e in Ei = e {aij{Bj) + 'r]ijk{Bj) ^j.). The factor e controls the weight of the turbulent EMF 
in relation to the differential rotation. 

• The ignorance of the a-effect in the polar regions is parameterized with with a ~ 
au cos 9 sin^ 6. 

• a exists only in the overshoot region, r/x only in the convection zone. 

• The correlation time is taken from Tcorr D/ut — 10^ s so that 17* ~ 5 results. 

• The tensors a and r/x are computed for a rms velocity profile by Stix (1991), the rotation 
law is directly taken from Christensen-Dalsgaard and Schou (1988). 

The main results from these models are: 

• The cycle period has just the same sensitivity to the BL's thickness as the shell dynamo, i.e. 
Tcyc[yr] ~ D/e for D in Mm {D ~ 15. . . 35 Mm). 

• Due to the rotational r/x-quenching the linear model yields the correct cycle time. Generally 
the 22-year cycle can only exist for a dilution of e ~ 0.5 (cf. Fig. ^). 

• For au — the magnetic activity is concentrated near the poles, for Ou — 1 it moves to the 
equator. 

• For too thin BL's there are too many toroidal magnetic belts in each hemisphere (Fig. ^. 
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Figure 8: Cycle period Tcyc as a function of the BL thickness D/R and and the dilution factor 
e for Ou = 0. LEFT: e = 0.5, the results are well represented by Tcyc = D [yr/Mm]. RIGHT: 
D = 35 Mm. The horizontal line indicates the solar cycle period of 22 years. 

The dependence of the cycle period on e is shown in Fig. ||. For e « 1 the linear solutions have 
the correct 22-year solar cycle period (right panel). The cycle period increases linearly with D, 



'^cyc 



(32) 



(see left panel of Fig. |^, which is in agreement with the linear relation (^ for spherical shell 
dynamos. For too thin boundary layers the cycle time will thus become very short. Nonlinearity 
further shortens the cycle period and this effect has to be compensated for by taking a smaller 
value of e, depending on whether or not magnetic buoyancy is included. 

In those cases where a-quenching and magnetic buoyancy are included, the rms-values of the 
poloidal and toroidal fields are 1 1 Gauss and 5 kGauss. Poloidal and toroidal fields as well as their 
mean and fluctuating parts are difficult to disentangle observationally. A toroidal field of 5 kGauss 
is, however, consistent with the observed total flux of 10^^ Mx, distributed over 50 Mm in depth 
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Figure 9: The magnetic field geometry for tfie nonlinear (a-quenclied) solution without buoyancy 
and e = 0.5, D = 35 Mm, = 1. An obvious problem of the BL dynamo is the large number of 
toroidal field belts. 

and 400 Mm in latitude. On the other hand, a poloidal field of 11 Gauss is larger than the observed 
value. 

In all cases the magnetic field exceeds the equipartition value which is here 6 kGauss. This is partly 
due to the turbulent diamagnetism leading to an accumulation of magnetic fields at the bottom of 
the convection zone, and partly due to too optimistic estimates for a-quenching. 
The dilution factor e is the main free parameter that must be chosen such that the 22-year magnetic 
cycle period is obtained. The thickness of the overshoot layer must be chosen to match the correct 
number of toroidal field belts of the Sun. Our results suggest that the thickness should not be much 
smaller than D ^ 35 Mm ^ 1 /2Hp. 

Nonlinear effects typically lead to a reduction of the magnetic cycle period. Unfortunately, there 
are no consistent quenching expressions that are valid for rapid rotation. Magnetic buoyancy also 
provides a nonlinear feedback (e.g. Durney and Robinson, 1982; Moss et al, 1990), but compu- 
tations by Kitchatinov and Pipin (1993) indicate that this only leads to a mean field advection of a 
few m/s, i.e. less than the advection due to the diamagnetic effect. However, magnetic buoyancy 
is important in the upper layers (where the diamagnetic effect is weak), and can have significant 
effects on the field geometry and the cycle period (Kitchatinov, 1993). 




4 Distributed Dynamos with Meridional Circulation 

Not only differential rotation but also meridional flow will influence the mean-field dynamo. 
This influence can be expected to be only a small modification if and only if its characteristic time 
scale Tdrift exceeds the cycle time Tcyc of about 11 years. With Tdrift — R/u^ we find ~ 2 
m/s as a critical value. If the flow is faster (Fig. 10) the modification can be drastic. In particular. 



if the flow counteracts the diffusion wave, i.e. if the drift at the bottom of the convection zone is 
polewards, the dynamo might come into trouble. The behaviour of such a solar-type dynamo (with 
superrotation at the bottom of the convection zone and perfect conduction beneath the convection 
zone) is demonstrated in the present Section. 



14 




Figure 10: The meridional circulation at mid-latitudes (45°) at the bottom (LEFT) and the top 
(RIGHT) of the convection zone for both a solar model and a K5 main-sequence star as a function 
of rotation period (Kitchatinov and Riidiger, 1999). The circle represents the solar case. The flow 
at the bottom is equatorwards (as it is the flow at the surface, there are 2 cells in radius). 



In the following the dynamo equations are given with the inclusion of the induction by merid- 
ional circulation. For axisymmetry the mean flow in spherical coordinates is given by u = 

{ur,U0,rsm6 fl). 

With a similar notation the magnetic field is 

/ 1 dA 1 dA \ 
V sm t) od r sin or ) 

with A as the poloidal-field potential and B as the toroidal field. Their evolution is described by 

dA , ^ „ (Q'^^ sine a / 1 

- + (u.V)A = -sB^r,^y-^^-^-[—--)y (34) 
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for uniform t^t and with s = rsm6 (cf. Choudhuri et al, 1995). In order to produce here a 
solar-type butterfly diagram with u = the a-effect must be taken as negative in the northern 
hemisphere (Steenbeck and Krause, 1969; Parker, 1987). Our rotation law and the radial merid- 



ional flow profile at 45° are given in Fig.[ll|. A one-cell flow pattern is adopted in each hemisphere. 
The approximates the latitudinal drift close to the convection zone bottom, the circulation is 
counterclockwise (in the first quadrant) for positive u™, i.e. towards the equator at the bottom of 
the convection zone. Note that the cycle period becomes shorter for clockwise flow (Roberts and 
Stix, 1972) and it becomes longer for the more realistic counterclockwise flow. The latter forms a 
poleward flow at the surface in agreement with observations and theory (Fig. [I^). For too strong 
meridional circulation {u™ > 50 m/s) our dynamo stops operation. The same happens already for 

— 10 m/s if the flow at the bottom is opposite to the magnetic drift wave. 
However, it is hard to explain why the convection zone itself is producing the negative a-values 



used in the model. More easy to understand are the positive values adopted in Fig. [13| (see Durney 
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Figure 11: Rotation law and the pattern of a meridional flow at 45° for the Steenbeck-Krause 
dynamo model of Section 4. The meridional circulation is prescribed (Riidiger, 1989) and not yet 
a result of differential rotation theory. The rotation profile always has a positive slope. 
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Figure 12: Operation of a Steenbeck-Krause dynamo with negative alpha-effect (in northern 
hemisphere) with and without meridional circulation. Eddy diffusivity is 5 • 10^^ cm^/s. Pos- 
itive bottom drift means counterclockwise flow (equatorwards at the bottom of the convection 
zone, polewards at the surface). TOP: The cycle period (solid, in years) and the magnetic energy 
(dashed) vs. the latitudinal drift at the bottom of the convection zone. MIDDLE: a-profile in ra- 
dius and butterfly diagram of the circulation-free dynamo. The maximal toroidal field is plotted vs. 
the time in years. BOTTOM: Butterfly diagrams for clockwise and counterclockwise circulations. 
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Figure 13: The same as in Fig. [12| but for positive a-effect (in northern hemisphere). Note the 
decay of the dynamo already for 10. . . 20 m/s (equatorwards at the bottom, polewards at the top of 
the convection zone). 



1995, 1996). Of course, now without circulation the butterfly diagram becomes 'wrong', with a 
poleward migration of the toroidal magnetic belts. One can also understand that the counterclock- 
wise flow (equatorwards at the convection zone bottom) rapidly destroys the dynamo operation, 
here already for 10 m/s amplitude. This amplitude is so small that it cannot be a surprise that such 
a dynamo might not be realized in the Sun. For the eddy diffusivity 5 • 10^^ cm^/s is used here, so 
the magnetic Reynolds number Rm = u^R/rn^ is 

Rm = 14-— — , (36) 
1(J m/s 

i.e. Rm = 14 for = 10 m/s. In their paper Choudhuri et al. (1995) apply Reynolds numbers of 
order 500 while Dikpati and Charbonneau (1999) even take Rm = 1400. Indeed, their dynamos 
are far beyond the critical values of Rm ~ 10 which are allowed for the weak-circulation aO- 



dynamo. Their dynamos exhibit the same period reduction for increasing shown in Fig. 13 but 
the butterfly diagram is opposite. It seems that quite another branch of dynamos would exist for 
Rm » 10. This is indeed true. 
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Figure 14: TOP: Butterfly diagram of a marginal dynamo with positive a-effect (3 m/s), with 
eddy diffusivity of 10^^ cm^/s and a bottom flow of 7 m/s. BOTTOM: {Br){B^), black colour 
indicates negative sign, i.e. the observed phase lay between poloidal and toroidal field. 



Table 1 : Cycle periods in years of solar-type-butterfly dynamos (bold) with positive a-effect (= 3 
m/s) and meridional flow of u™ at the bottom of the convection zone, the eddy diffusivity given in 
the headline. Signature ' — ' indicates polewards migration, 'decay' means decaying fields. Note 
the good representation of the r/-dependence of the cycle time of @ in the first line of the Table. 
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The butterfly diagram of such a dynamo with positive a-effect (3 m/s) and a magnetic Reynolds 
number of 450 is given in Fig. 14. It works with low eddy diffusivity of 10^^ cm^/s and the flow 



pattern of Fig. |11| In the lower panel the phase relation of poloidal and toroidal magnetic fields 
is shown. In lower latitudes the observed phase lay between the field components is obtained. In 
Table D the cycle periods in years are given for positive-a dynamos. Bold numbers indicate the 
solar-type butterfly diagram. For fast meridional flow and small eddy diffusivity we indeed find 
oscillating dynamos. The eddy diffusivity is the same as that known from the sunspot decay, i.e. 



10^^ cm^/s (see Brandenburg, 1993). 
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5 The Maunder Minimum 



The period of the solar cycle and its amplitude are far from constant. The most prominent activity 
drop was the Maunder minimum between 1670 and 1715 (Sporer, 1887). The variability of the 
cycle period can be expressed by the 'quality', cJcyc/Awcyo which is as low as 5 for the Sun 
(Hoyng, 1993, cf. Fig. |). 

Measurements of ^^C abundances in sediments and long-lived trees provide much longer time 
series than sunspot counts. In agreement with Schwarz (1994), Vos et al. (1996, 1997) found a 
secular periodicity of 80-90 years as well as a long-duration period of about 210 years. The mea- 
surements of atmospheric ^^C abundances by Hood and Jirikowic (1990) suggested a periodicity 
of 2400 years which is also associated with a long-term variation of solar activity. The variety of 
frequencies found in solar activity may even indicate chaotic behaviour as discussed by Rozelot 
(1995), Kurths et al. (1993, 1997) and Knobloch and Landsberg (1996). 

The activity cycle of the Sun is not exceptional: The observation of chromospheric Ca-emission 
of solar-type stars yields activity periods between 3 and 20 years (Noyes et al., 1984; Baliunas 
and Vaughan, 1985; Saar and Baliunas, 1992a,b). A few of these stars do not show any significant 
activity. This suggests that even the existence of the grand minima is a typical property of cool 
main-sequence stars like the Sun. From ROSAT X-ray data Hempelmann et al. (1996) find that 
up to 70% of the stars with a constant level of activity exhibit a rather low level of coronal X-ray 
emission. HD 142373 with its X-ray luminosity of only log Fx = 3.8 is a typical candidate. We 
conclude that during a grand minimum not only the magnetic field in the activity belts is weaker 
than usual but the total magnetic field energy is also reduced. 

The short-term cycle period appears to decrease at the end of a grand minimum according to a 
wavelet analysis of sunspot data by Frick et al. (1997). There is even empirical evidence for a 
very weak but persistent cycle throughout the solar Maunder minimum as found by Wittmann 
(1978) and recently Beer et al. (1998). The latitudinal distribution of the few sunspots observed 
during the Maunder minimum was highly asymmetric (Sporer, 1887; Ribes and Nesme-Ribes, 
1993; Nesme-Ribes et al, 1994). Short-term deviations from the north-south symmetry in regular 
solar activity are readily observable (Verma, 1993), yet a 30-year period of asymmetry in sunspot 
positions as seen during the Maunder minimum remains a unique property of grand minima and 
should be associated with a parity change of the internal magnetic fields. 

5.1 Mean- field Magnetohydrodynamics 

The explanation of grand minima in the magnetic activity cycle by a dynamo action has been 
approached by two concepts. The first one considers the stochastic character of the turbulence 
and studies its consequences for the variations of the a-effect and all related phenomena with 
time (see Section ||). The alternative concept includes the magnetic feedback to the internal solar 
rotation (Weiss et al., 1984; Jennings and Weiss, 1991). Kitchatinov et al. (1994a) and Tobias 
(1996, 1997) even introduced the conservation law of angular momentum in the convection zone 
including magnetic feedback in order to simulate the intermittency of the dynamo cycle. 
A theory of differential rotation based on the A-effect concept is coupled with the induction equa- 
tion in a spherical 2D mean-field model. The mean-field equations for the convection zone include 
the effects of diffusion, a-effect, toroidal field production by differential rotation and the Lorentz 
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force. They are 



dA d'^A sine d f 1 dA\ 



dB _ Id/ d{Br)\ tjt d / 1 d{Bsme)\ 
dt r dr \^ dr J r"^ d6\ sin 6 dO ) 

ld9.dA IdVidA 1 d f dA\ I d / a dA\ 



r dr 86 r dO dr rs'mO dr \ dr J r^ dO \sm9 dO J ' 

^ 1 ndAd{Br) 1 dAd{B sin e)\ ^^^^ 
fiQr"^ sin 6 \r do dr sin 6 dr dO J 

Since the importance of the large-scale Lorentz force in the momentum equation was discussed by 
MaUcus and Proctor (1975), we call this term 'Malkus-Proctor' effect. 

The computational domain is a spherical shell covering the outer parts of the Sun down to a 
fractional solar radius, x = r/R, of 0.5. The convection zone extends from x = 0.7 to x = 1. 
The a-effect works only in the lower part from x = 0.7 to x = 0.8 while turbulent diffusion 
of the magnetic field, turbulent viscosity, and the A-effect are present in the entire convection 
zone. Below x = 0.7 both the magnetic diffusivity and the viscosity are two orders of magnitude 
smaller than in the convection zone. The boundary conditions are specified as Qr^ = dA /dr = 
B = at r = R, and Qr,j, = A = B = at the inner boundary. The angular momentum 
flux is given by 

Qr^ = VT sine (^-r ^ + , Qe^ = -v^ sin^ ^. (40) 

y^^'^ determines the radial rotation law without magnetic field. 

Five dimensionless numbers define the model, namely the magnetic Reynolds numbers of both 
the differential rotation and the a-effect, Cq, = Q.qR^ /rjT, Ca = aQR/r]T, the magnetic Prandtl 
number Pm = i^t/^?t ; the Elsasser number 

B^ 

E = V; (41) 

and the A-effect amplitude, V^^h In the a-effect, 

a = ao cos 9 sin^ 6, (42) 

the factor sin^ 6 has been introduced to restrict magnetic activity to low latitudes and ao — 
Icovv ^o> SO that Cfi/lCal — R/lcorv, hcncc in general Cq exceeds Ca ('aO, dynamo'). Our 
dynamo works with = — 10 and Cq = 10^. V^'^^ is positive in order to produce the required 
super-rotation, and its amplitude is 0.37. With the eddy diffusivity (|1|) the Elsasser number reads 
E = IjcrfL* and is set to unity here. (See also Kiiker et al, 1999; Pipin, 1999.) 
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Figure 15: LEFT: The time dependence of the dynamo for the large-scale Lorentz-force feedback 
only (Pm = 0.1, E = 1). Top: Toroidal magnetic field (left) and magnetic energy (right). Bottom: 
Cycle time (left) and magnetic parity (right). RIGHT: The same with a-quenching but still without 
A-quenching (A = 0). 
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Figure 16: The same as in Fig. |T^ but for strong A-quenching (A = 25). TOP: Toroidal magnetic 
field (left) and magnetic energy (right). BOTTOM: Cycle time (left) and magnetic parity (right). 
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5.2 Results of the Interplay of Dynamo and Angular Velocity 



Figs. |15| and |16| demonstrate the action of different effects and show the variation of the toroidal 
magnetic field at a fixed point (x = 0.75, 9 = 30°), the total magnetic energy, the variation of the 
cycle period, and the parity P = {Es — Ea)/{Es + Ea), derived from the decomposition of the 
magnetic energy into symmetric and antisymmetric components (Brandenburg et ai, 1989). All 
times and periods are given in units of a diffusion time E? /rn^. Field strengths are measured in 
units of i?eq- 

If the large-scale Lorentz force ('Malkus-Proctor effect') is the only feedback on rotation, the time 
series given in Fig. |^ may be compared with the results in Tobias (1996) although a number of 
assumptions are different between Tobias' Cartesian approach and our spherical model. Fig. p3| 
shows a quasi-periodic behaviour with activity interruptions like grand minima. This model, how- 
ever, neglects the feedback of strong magnetic fields on the a-effect and the differential rotation. 



The right-hand side of Fig. 15 shows the result of the same model but with a local a-quenching. 



where i^tot is the absolute value of the magnetic field. The variability of the cycles turns into a 
solution with only one period. Similar to the suppression of dynamo action, a quenching of the 
A-effect causing the differential rotation is according to 

If A is near unity, the maximum field strength and total magnetic energy decrease slightly, but the 
periodic behaviour remains the same, i.e. the effect of the A-quenching is too small to alter the 
differential rotation significantly. However, an increase of A leads to grand minima - an example 
for A = 25 is given in Fig. [l^. Minima in cycle period occur shortly after a grand activity 
minimum in agreement with the analysis of sunspot data by Frick et al. (1997a,b). The amplitude 
of the period fluctuations in Fig. [l^ is much lower than in the Malkus-Proctor model but is still 
stronger than that observed. The magnetic Prandtl number used for the solutions in Figs. |l^ and 



Iq was Pm = 0.1 and it is noteworthy that grand minima do not appear for Pm = 1. 
The strong variations of the parity between symmetric and antisymmetric states appearing in the 
nonperiodic solutions can be explained by the numerous (usually 5 ... 6) magnetic field belts mi- 
grating towards the equator. Slight shifts of this belt-structure against the equator result in strong 
variations in the parity. Averaged over time, dipolar and quadrupolar components of the fields 
have roughly the same strength; all periodic solutions have strict dipolar structure. 



Spectra of long time series of the toroidal magnetic field are given in Fig. |17| for both the Malkus- 
Proctor model and the strong A-quenching model. The long-term variations of the field will be 
represented by a set of close frequencies whose difference is the frequency of the grand minima. 
The Malkus-Proctor model shows a number of lines close to the main cycle frequency. The dif- 
ference between the two highest peaks can be interpreted as the occurrence rate of grand minima. 
However, the shape of the spectrum indicates that the magnetic field appears rather irregularly. The 



spectrum of the model with all feedback terms and strong A-quenching is also given in Fig. 17 and 
shows a similar behaviour with highest amplitudes near the main cycle frequency of the magnetic 
field. The average frequency of the grand minima is represented by the distance between the two 
highest peaks. 
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Figure 17: LEFT: Power spectrum of the magnetic-field amplitude variations for the Malkus- 



Proctor model of Fig. 15. The frequency is given in arbitrary units. RIGHT: The same but for 
the model with strong A-quenching. The highest peaks are the main cycle frequency, whereas the 
difference between two peak neighbours indicates the secular cycle. 




Figure 18: LEFT: Correlation between magnetic field oscillations and variations of the differential 
rotation measure (dQ/dr)^, averaged over the latitude 6. RIGHT: Correlations between magnetic 
field energies for both parities (dipolar and quadrupolar energy component) and the differential 
rotation. For most of the time the orbit resides in the state of regular cycles in the upper part of the 
box. Differential rotation is suppressed during grand minima. 
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The inteiplay of magnetic fields and differential rotation is demonstrated in Fig. 18 showing the 
variation of the radial rotational shear averaged over latitude, {{dU, / dr)"^) , versus time, compared 
with the toroidal field. A minimum in differential rotation is accompanied by a decay of the 
magnetic field and followed by a grand minimum. The differential rotation is being restored 
during the grand minimum since the suppressing effect of magnetic A-quenching is reduced. The 
behaviour of the system can be evaluated in a phase diagram of the dipolar component of the 
magnetic field energy, the quadrupolar component £'s, and the mean angular velocity gradient, 
as also shown in Fig. 18. The dots in the interior of the diagram represent the actual time series; 
projections of the trajectory are given at the sides. 

The trajectory resides at strong differential rotation during normal cyclic activity. The magnetic 
field oscillates in a wide range of energies. The trajectory moves down as the field starts to suppress 
the differential rotation, with oscillation amplitudes diminishing. The actual grand minimum is 
expressed by a loop at very low dipolar energies; the quadrupolar component, however, remains 
present throughout the minimum. The differential-rotation measure is already growing at that 
time. At first glance, the phase graph may be associated with Type 1 modulation as classified in 
Knobloch et al. (1998). 

The effect of large-scale Lorentz forces on the differential rotation as the only feedback of strong 
magnetic fields produces irregular grand minima with strong variations in cycle period. The 
complex time series turns into a single-period solution, if the suppression of dynamo action (a- 
quenching) is included. If a strong feedback of small-scale flows on the generation of Reynolds 
stress (A-quenching) is added, grand minima occur at a reasonable rate between 10 and 20 cycle 
times. The cycle period varies by a factor of 3 or 4. The northern and southern hemispheres dif- 
fer slightly in their temporal behaviour. This is a general characteristic of mixed-mode dynamo 
explanations of grand minima. 
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The magnetic Prandtl number directs the intermittency of the activity cycle. Values smaller than 
unity are required for the existence of grand minima, whereas the occurrence of grand minima 
again becomes more and more exceptional for very small values of Pm (Fig. 19). 



6 Stellar Cycles 

The cycle period may be considered as a main test of stellar dynamo theory because it reflects 
an essential property of the dynamo mechanism. A realistic solar dynamo model should provide 
the correct 22-year cycle period and, in the case of stellar cycles, the observed dependence of the 
cycle period on the rotation period. 

The following Sections will deal with a number of applications of ID models illuminating some 
effects of more complicated EMF. This will include the additional consideration of the magnetic- 
diffusivity quenching, a dynamo-induced a-function, and the influence of temporal fluctuations of 
a and r^T- 

The relations between stellar parameters and the amplitude and duration of magnetic activity cy- 
cles are fundamental for our knowledge of stellar physics. Observational results cover a variety of 
findings which depend on how stars are grouped and how unknown parameters must be chosen. 
Generally, relations for the magnetic field and the cycle frequency are expressed by 

S^ax OC O*"*, Weyc OC (45) 

with different m and n. The value discussed by Noyes et al. (1984) was n = 1.28. Saar (1996) 
derives a very small m (his Fig. 3), while Baliunas et al. (1996) find n ~ 0.47 for young stars and 
n ~ 1.97 for old stars. Brandenburg et al. (1998) derive nj = 1.46 and toa = 1-48 for inactive 
(I) and active (A) stars. Saar and Brandenburg (1999) find n = 0.6 for the young ("super-active") 
stars and n = 1.5 for the old ones. Higher values are reported by Ossendrijver (1997) consistent 



with n = 2 . . . 2.5. All the reported exponents so far are positive. The relations in (45) can be 
reformulated with the dynamo number V as 

Smax = BeqP™', Wcyc = ^^^'/^diff, (46) 

where i?cq means the 'equipartition field' and Tdiff = H"^ /rjQ the diffusion time. Instead of the 
normalization used in (p6|)2. Soon et al. (1993) proposed the use of the basic rotation rate a 
procedure we do not follow here to retain consistency with conventional definitions of the dynamo 
number. 

The main issues of stellar activity and the relations to dynamo theory concerning cycle times were 
formulated by Noyes et al. (1984). Their argumentation concerns a zeroth-order aO-dynamo 
model for which the coefficient n in the scaling cjcyc oc D"/^ is determined. The linear case yields 
n = 4/3 for the most unstable mode (Tuominen et al., 1988). Different nonlinearities lead to 
different exponents. For a-quenching, n = is found while for models with flux loss (rather 
than a-quenching) n = 1 results or n = 2/3 - if only the toroidal field is quenched by magnetic 
buoyancy. Note the absence of finite n for simple a-quenching. If the theory is correct and the 



observations are following relations such as (45), then the nonlinear dynamo can never work with 
a-quenching as the basic nonlinearity. 

A number of more complex dynamo models confirming the findings of Noyes et al. (1984) have 
been developed in the last decade. For large dynamo numbers Moss et al. (1990) also derived n = 
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1 from a spherical dynamo saturated by magnetic buoyancy. Schmitt and Schiissler (1989) showed 
with a ID model (their Fig. 6) that there is practically no dependence of the cycle frequency on the 
(large) dynamo number for a-quenching, but they find a strong scaling (n ~ 2) for their 'flux-loss 
models'. Also for a 2D model in spherical symmetry the dependence of the cycle frequency on 
the dynamo number proved to be extremely weak, n < 0.1 (Riidiger et al, 199 A, dashed line in 
their Fig. 4). 

The latter application, however, requires further consideration. One finds a finite value for n 
if the magnetic feedback is not only considered acting upon the a-effect, but also for the eddy 
diffusivity tensor. What we assume here is that the magnetic field always suppresses and deforms 
the turbulence field, and this has consequences for both the a-effect and the eddy diffusivity. 
The theory for a second-order correlation-approximation is given in Kitchatinov et al. (1994b), 
applications are summarized in Riidiger et al. (1994). 

The turbulent-diffusivity quenching concept was also described in Noyes et al. (1984). As the 
magnetic fields become super-equipartitioned, however, a series expansion such as used cannot 
be adopted. Tobias (1998), with a 2D global model in Cartesian coordinates, finds n' varying 
between 0.38 and 0.67 depending on the nonlinearity. He also finds the exponent growing with 
increasing effect of j^-quenching. The results for the Malkus-Proctor effect alone yield the weakest 
dependence for the cycle period with V. 

Since the paper by Brandenburg et al. (1998), 'anti-quenching' expressions such as 

a oc BP, r/T oc 5« (47) 

are also under consideration, with p > and q > {). The idea is to find the consequences of EMF 
quantities being induced by the magnetic field itself (see Section 7). Saturation of such dynamos 
requires q > p. All their global models, however, lead to exponents n of order 0.5. It concerns 
both ID and 2D models, characteristic values for q were numbers up to 6. 

It is shown in the following that a special ID slab dynamo, as a generalization of Parker's zero- 
dimensional wave dynamo without r^-quenching, has a very low n while n exceeds unity if -q- 
quenching is included. Buoyancy effects are therefore not the only possible addition to dynamo 



theory to explain observations as given in (45 ); the magnetic suppression of the turbulent magnetic 



diffusivity also yields an explanation of the observations. 
6.1 The Model 

Here we are not working with the simplest possibility, which would read a-ij = adij and rjij^ = 
VT^ijk- The full feedback of the induced magnetic field on the turbulent EMF is included, i.e. the 
magnetic suppression and deformation of both the tensors a and r]. In particular, the influence 
of the magnetic field on the magnetic diffusivity is often ignored in dynamo computations, and 
we shall demonstrate the differences in the gross properties of the solutions with and without rj- 
quenching. The main consequence of the inclusion of r/-quenching is the appearance of a nonlinear 
'magnetic velocity' in 6, 

£ = ■■■ + U™^s X (B) (48) 

with 

U- = OVlog(i?) + ,.^^^^i?I^ (49) 
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Figure 20: The cycle frequency of oscillatory solutions of the a^il^-dynamo versus Cq with 
different a-quenching functions. A second solution can be excited for the slower cut-off func- 
tion (dashed) for high Cq. The steeper cut-off function (sohd) delivers chaotic solutions above 
Cn ~ 200. 
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The components a,^^ and ass of the a-tensor are taken from Riidiger and Schultz (1997). They 
are equal in our approximation and read 



(55) 



with the a-quenching function ^' given in (|29|). These quenching functions are approximated by 
*(/3) ~ 1 - fp^, (p{(3) ~ 1 - f /32, ~ I/32 and ^(/3) ~ f /32 for weak magnetic fields. 

The only coordinate about which our quantities are allowed to vary is the direction of the rotation 
axis, i.e. z. The z-dependencies in (55) are summarized in the form a = —a sin2z where the 
lower and upper boundary are located at z = and z = ir. The plane has infinite extent in the 
X and y directions and is restricted by boundaries in the z-direction. The normalized dynamo 
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Figure 21: LEFT: The cycle frequency of oscillatory solutions of the a^Q-dynamo vs. Cq with 
(solid) and without (dashed) r/-quenching. RIGHT: The exponent n' for the scaling (36). Only 
diffusivity quenching yields exponents of order 0.5. The model without r/-quenching contains a 
second solution which can be excited at high Cq. 



equations are 

dA ^ ,rr. dA 

dB ^ d f , , dA\ ^ dA d f dB\ d / n 

with w = {2(p — ipz)dlog Btot/dz and B'^^^ = B"^ + {dA/ dz)"^ . The boundary conditions are i? = 
dA /9z = at z = and z = vr which limits the magnetic field to the disk. The dimensionless 
parameters 

a = Cn = (57) 

represent the strength and sign of the a-effect and the shear flow. Ca is fixed to a value of 5 here. 
A positive Ca describes a positive a in the upper half of the layer and vice versa. We chose to 
vary Cq with positive values which represent a positive shear. 



6.2 The Results 

All the magnetic fields of the nonlinear ID model are symmetric with respect to the equatorial 
plane. In our domain of positive shear the solutions of both the a^-dynamo (small Cq) and the 
aO-dynamo (large Cq) are oscillatory. The transition between these regimes occurs at Cq — 50 



(cf. Riidiger and Arlt, 1996). The cycle periods are always constant over time. Fig. gO| shows 
the resulting cycle frequencies for the conventional a-quenching expression. The cycle period is 
strikingly independent of the dynamo number (see also Noyes et ai, 1984; Schmitt and Schiissler, 
1989; Jennings and Weiss, 1991). On the other hand, Fig. ^gives the resulting cycle frequencies 
of the oscillatory solutions for the a^il-dynamo with the /^-quenching concept. The cycle periods 
are given for models with and without r/-quenching. The exponent n' is variable and positive for 
the model with ry-quenching whereas it is much more complex for the model without 77-quenching. 
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The resulting exponent n' ~ 0.5 for large dynamo numbers in Fig. 21 is remarkably similar to 
the behaviour of linear dynamos. In order to compare this value with the observed value of n in 
(4^) we have to introduce the scaling of the physical quantities with il*. We adopt the scaling 
r/o oc from Kitchatinov et al. (1994b) and use the parameterization 



d\ogr 



oc f)*". (58) 



Negative k describe a decrease of the unknown normalized differential rotation in stellar convec- 
tion zones with increasing basic rotation. Negative k are indeed produced by theoretical models 
of the differential rotation, though only if meridional flow is included. The results in Fig. ^lead to 
K ~ — 1. The dynamo number scales as 

V oc (59) 

if the a-effect is assumed as scaling with W^^. Then the relations (45) and (p^) yield 

n = + 3 + K)n - 1. (60) 

We apply the exponent n' ~ 0.5 from Fig. |1| for large dynamo numbers and get n = (/i + k)/2 + 
0.5, so that with k = — 1 we find ^ ~ 2n. Consequently, the observed exponents n = 0.5 ... 1.5 
would lead to 

^=1...3 (61) 
for the a-effect running with the Coriolis number which does not seem too unreasonable. 



7 Dynamo-induced on-off Alpha-Effect 

A two-component model for the solar dynamo has been suggested by Ferriz-Mas et al. (1994). 
Toroidal magnetic fields are generated by differential rotation in the overshoot layer whereas a 
dynamo acts in the upper part of the convection zone. The thin overshoot layer is assumed to be 
stably stratified up to magnetic field strengths of 10^ Gauss (Ferriz-Mas and Schiissler, 1995). The 
dynamo may operate in the turbulent convective zone and impose stochastic perturbations to the 
overshoot layer. 

Schmitt et al. (1996) discuss whether the interaction of the weak-field convection zone dynamo 
and the magnetic fields in the overshoot layer might explain the grand-minimum behaviour of 
the solar cycle. They assumed a mean-field dynamo with an a-effect working only when the 
magnetic field exceeds a certain threshold. Grand minima are expected when the perturbations 
from the convection zone do not keep the overshoot layer instabilities above this threshold. Here 
we test whether a dynamo, once excited, will really turn off in this sense. 

The dynamo equations are studied in both ID and 2D domains. In the ID approach the integration 
region extends only in the z-direction whence the remaining radial and azimuthal components of 
the magnetic field depend on z only. Normalization of times with the diffusion time Tdiff = H'^ /rj^r 
and vertical distances z with H yields 

8 A 8"^ A 

— = CMz)^iBtot)B + ^, (62) 
dB ^ d f.,.^,^ .dA\ ^ dA d^B 
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Figure 22: The toroidal component of the magnetic field (LEFT) at a fixed point of the domain 
for two different on-off a-functions (RIGHT) in a ID model. The initial conditions imply a field 
strong enough to cover partly the 'on'-range of a. 
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Figure 23: The toroidal component of the magnetic field (LEFT) for two different on-off a- 
functions (RIGHT) in a 2D model. Smax = 1 in both cases. 
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The a depends on the magnetic field as well as on the location in the object. Hence we have two 
different dependencies in a = a(z)ao^(-Btot)- Again, the first factor is expressed by a{z) = 
— sin 2z. 

The a-quenching by magnetic fields has usually been expressed by functions such as 1/[1 + 
(-Btot/^max)^]^ with a cut-off field strength -Bmax related to the energy of velocity fluctuations 
or the gas pressure. In the present approach a acts only in a certain range of i?tot and it is zero 
otherwise, i.e. 

The lower threshold -Bmin represents the onset of instability in the strong-field dynamo layer at 
the bottom of the convection zone, and the upper threshold -Bmax is related to the buoyancy-driven 
disappearance of flux tubes from the overshoot layer. Again a vacuum surrounds the ID model. 
The normalized induction equations of the 2D model read 

dB _ Id/ d{Br)\ t]^ d / 1 d{Bsme)\ 
dt r dr dr ) r^d9\ sin 6 89 J 

^ IdndA 1 d f^dA\ I d / a dA\ ^^^^ 
r dr 86 rsm.9 Or \ dr J 09 \sin9 89 J ' 

OA 8"^ A sm9 8 ( 1 8A\ 

-Q^ = ^^ + ^^^00 [-^9m) +«"-^^- (66) 

The computational domain extends from r = 0.5 to r = 1.0 measured in solar radii, and from 
9 = to 9 = TT in colatitude. The dynamo effect is assumed to be placed between r = 0.7 and 
r = 0.8 neighboured by a weakly dissipative layer with t/t = 0.01 below and a zone with strong 
turbulent diffusion (?7t = 1) above the dynamo layer; rj^ is also unity in the dynamo shell. Here 
we neglect the ^-dependence of the angular velocity and assume a constant 80. / 8r for simplicity. 
An a^il-dynamo provides oscillatory and steady-state solutions. The first test in a ID model 
applies Cq = 100 and Cq = 5 and two on-off functions with both i?min = 0.5 and -Bmin = 0.9. 
The dynamo was excited with an initial magnetic field strong enough to reach the 'on' -range of a. 



Fig. 22 shows the resulting oscillatory solutions. The dynamo never dies since the a-effect acts 
before the magnetic fields completely 'dives' through the on-off function. 

The same is found in a 2D model as shown in Fig. We used the parameters Cq = 10^ and 
Ca = —10 and two on-off functions with i?min = 0.38 and i?min = 0.5, resp. If the on-off 
function is too narrow the solutions change from oscillatory to steady-state, although we never 
find the dynamo stopping its operation. The parity of the oscillatory solution is antisymmetric 
with respect to the equator but equator symmetry is found for the steady-state solution. 
The results indicate that the dynamo, once initiated by sufficiently strong magnetic fields, will 
not switch off on its own and, hence, does not impose grand minima to its cyclic behaviour. The 
dynamo will only cease operating if the flux tubes migrate into the upper convection zone on a 
time scale much shorter than the magnetic diffusion time. 



8 Random Alpha 

Two different types of magnetic dynamos appear to be acting in the Sun and the Earth. In contrast 
to the distinct oscillatory solar magnetic activity, the terrestrial magnetic field is 'permanent'. 
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Indeed, two different sorts of dynamos in spherical geometry have been constructed with such 
properties. The a^-dynamo yields stationary solutions while the aO-dynamo yields oscillatory 
solutions. Observation and theory seem to be rather in agreement. There are, however, differences 
in the time behaviour: solar activity is not strictly periodic and the Earth's dynamo is not strictly 
permanent. The Earth's magnetic field also varies. There are irregular changes of the field strength, 
the shortest period of stability being 40 000 years. The average length of the intervals of constant 
field is almost 10 times larger than this minimum value (cf. Krause and Schmidt, 1988). Such a 
reversal was simulated numerically with a 3D MHD code by Glatzmaier and Roberts (1995). 
In Section 5 we explained grand minima such as the Maunder minimum by various nonlinearities 
in the mean-field equations. The averages aie taken over an 'ensemble', i.e. a great number of 
identical examples. The other possibility to explain the temporal irregularities is to consider the 
characteristic turbulence values as a time series. The idea is that the averaging procedure concerns 
only a periodic spatial coordinate, e.g. the azimuthal angle In other words, when expanding in 
Fourier series such as e*"*"^ the mode m = is considered as the mean value. Again, if the time 
scale of this mode does not vary significantly during the correlation time, local formulations such 
as equation ( |69| ) below are reasonable. Nevertheless, the turbulence intensity, the a-effect, and the 
eddy diffusivity become time-dependent quantities (Hoyng, 1988; Choudhuri, 1992; Hoyng, 1993; 
Moss et at, 1992; Hoyng et ai, 1994; Vishniac and Brandenburg, 1997; Otmianowska-Mazur et 
al, 1997). 

However, questions concerning the amplitude, the time scales, and the phase relation between (say) 
helicity and eddy diffusivity need to be studied. Are the fluctuations strong enough to influence the 
dynamo significantly? We construct here a random turbulence model and evaluate the complete 
turbulent electromotive force as a time series. The consequences of this concept are then computed 
for a simple plane-layer dynamo with differential rotation. The main parameter which we vary is 
the number of cells in the unit length. 

One can define mean values of a field F by integration over (say) space. Of particular interest here 
is an averaging procedure over longitude, i.e. 



(Braginsky, 1964; Hoyng, 1993). The turbulent EMF for a given position in the meridional plane 
forms a time series with the correlation time Tcorr as a characteristic scale. The peak-to-peak 
variations in the time series should depend on the number of cells. They remain finite if the 
number of cells is restricted as it is in reality. For an infinite number of the turbulence cells the 
peak-to-peak variation in the time series goes to zero but it will grow for a finite number of cells 
along the unit length. 

The tensors constituting the local mean-field EMF in (^, a and rj, must be calculated from one 
and the same turbulence field. We propose to define a helical turbulence existing in a rectangular 
parcel and to compute simultaneously their related components. 

We restrict ourselves to the computation of the turbulent EMF in the high-conductivity limit. Then 
the second-order correlation approximation yields 




(67) 







oo 




(68) 
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which for short correlation times can be written in the form (Q). A Cartesian frame is used where 
y represents the azimuthal direction over which the average is performed. For the dynamo only 
the components and £y are relevant. 

It would be tempting to apply ( |68l ) as it stands. In components it reads 

d{By) ^ 
dz ' 



£x = Oixx{Bx)+rjT 



£y = ayy{By) -r]T^^^. (69) 



r/T plays the role of a common eddy diffusivity. From (|68|) one can read 



a..r = / {u'y{t) ) dr, 





oo 



riT = I {u',{t) u',{t - t)) dr (70) 
(Krause and Radler, 1980). 



8.1 The Turbulence Model 



We study the time evolution of the dynamo coefficients (|70|) generated by turbulent gas motions. 
We analyse a parcel of the solar gas situated in the convective zone. The rectangular coordinate 
system has the xy-plane parallel to the solar equator, the x-axis pointing parallel to the solar radius, 
the y-axis tangential to the azimuthal direction and the z-axis directed to the south pole. The parcel 
is permanently perturbed by vortices of the form of rotating columns of gas randomly distributed 
at all angles in 3D space. At every time step the resulting velocity field is used to calculate the 



EMF-coefficients after (70). 




Figure 24: Cross-section of the velocity field for the 100-eddy turbulence model B in the xy-plane 
at an arbitrary time of the simulation. White colour denotes positive Uz, black colour negative Uz- 
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A single vortex rotates with the velocity u around its axis of rotation and moves with the velocity 
up and down along this axis. The eddy itself has also its own lifetime growing Unearly with 
subsequent time steps. The vortex velocity uj as well as the vertical drift w decrease with the 
distance r from the axis and with the lifetime t according to 

-0.b\{r/lcorrf + {z/Zcorrf] ,. 

[io.w) = {loq.wo) ■ e V /.e-*/-^™", (71) 

where /corr is the adopted vortex radius (Otmianowska-Mazur et al, 1992), Zcorr is its length scale 
along the axis of rotation and Tcorr is the eddy decay time scale. All velocities are truncated at 
SZcorr and Szcorr- In the model all vortices are orientated Uke right-handed screws. The resulting 
motion has maximum helicity. A uniform distribution of eddies in a 3D space is approximated 
with 12 possible inchnation angles of their rotational axes to the a;y-plane. 
The initial state is a number TV of moving turbulent cells having random inclinations and positions 
in the xy-plane. After an assumed period of time (one or more time steps), a fraction of them, i2tur 
(the ratio of the new-eddies number to N), is changed to new ones with different randomly given 
position, inclination and with lifetime starting from zero. The situation repeats itself continuously 
with time. 

The vortex radius Zcorr as well as the decay time Tcorr and the fraction of the new turbulent cells 
i?tur are varied for four cases given in Tables 2 and 3 in normalized units (cgs units for time, 
velocity and diffusivity result after multiphcation with 2.5 • 10^ s, 10^ cm/s and 2.5 • 10^"^ cm^/s). 
We made numerical experiments with the following values for the vortex radius /corr and the length 
scale 2;corr in the 2;-direction in normalized units: 1, 2, 4 and 8. For successive length scales we 
choose the time scales 10, 20, 40 and 80. The ratio of both scales is always /corr /''"corr — 0.1. In 
order to obtain the expected mean turbulent velocity value of 0.1, the fraction of new vortices i?tur 
is also varied (Table 2; cf. Otmianowska-Mazur etal, 1997). 

8.2 Numerical Experiments 

The simulations deUver time series of the turbulence intensity, the eddy diffusivity and the a- 
coefficients, and a standard deviation a from their temporal averages 

a = E[{X - E{X))% (72) 

where E{X) is the time average of a random variable X. The standard deviation measures the 
amplitude of the fluctuations compared with the temporal mean of a given quantity. Table 2 lists 



Table 2: Input and output for the turbulence models v4, B, C and V. N is the eddy population of 
the equator, i?tur the eddy birth rate, other quantities and normalizations as explained in the text. 





/corr 


"^corr 


N 


-Rtur 


r/T 


o-(?7t) 




(y{.Olxx) 


ayy 


a{ayy) 


A 


1 


10 


200 


1.0 


0.191 


0.086 


-0.046 


0.048 


-0.011 


0.055 


B 


2 


20 


100 


0.7 


0.336 


0.208 


-0.029 


0.063 


-0.048 


0.059 


C 


4 


40 


50 


0.005 


1.045 


0.860 


-0.029 


0.119 


-0.039 


0.110 


V 


8 


80 


25 


0.0003 


1.793 


2.278 


-0.015 


0.136 


-0.017 


0.110 
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Table 3: The ratio S of fluctuations and mean values. 



model 




^Ctxx 




A 


0.45 


1.03 


0.72 


B 


0.62 


2.13 


1.24 


C 


0.82 


4.10 


2.79 


V 


1.27 


8.94 


6.62 



the input model parameters /corr, Tcorr, and Rtm and the coefficients obtained. The ratios of the 
standard deviations to their mean values, 

Sr, = air])/ \7]\, (73) 

are given in Table 3. 

The results for a model A with /corr = 1 and Tcorr = 10, which is the sample of shortest life- 
times and smallest spatial dimensions of individual vortices in our simulation sequence is given 



in Fig. |25|. All computed quantities are drawn as a time series, each in its own normalized units 
according to their definitions. The three EMF-coefficients as well as the turbulence intensity show 
fluctuations around their mean values in time. The diffusion coefficient r/x possesses positive 
values during most of the time. There are also short periods with a negative 7]^. The negative 
values are not significant because the resulting mean is always positive yielding 0.191 (Ta- 
ble 2) in normalized units. The ratio 5,,^ is 0.45 (Table 3). Indeed, the value 0.191 corresponds 
to r/T — Tcorr u^, which givcs 0.1 for this case. In contrast to ?7t, the coefficients Uxx and Uyy 
possess negative values during the major part of the time; for short periods, however, positive 
values are also present. The negative a results from the assumed non-zero, right-handed helicity 
of the vortices and is in agreement with the expectations. Averaged in time axx for model A is 
—4.6 • lO^cm/s. For ayy we get —7.7 • lO^cm/s resulting in a dynamo number of Ca = 10. It 
can be seen that axx is always slightly smaller than Uyy. This fact is connected with the assumed 
averaging only along the y-axis, which certainly influences the value of axx- Model A yields the 
ratios a{axx)/ \ axx \= 1-03 and a{ayy)/ \ Oyy \= 0.72. 

The second experiment, B, uses doubled correlation lengths and times, /corr = 2, Tcorr = 20 



(see also Fig. g4|). A fraction of 0.7 new vortices at every time step was applied. The resulting 
fluctuations of rjT and a-coefficients are higher than in the case A (Table 2). The averaged value 
of r/T increases to 0.336, yielding 8.5 • 10^^ cm^/s. The ratios Sa^^ and Sayy increase up to 2.13 
and 1.24, respectively (Table 3). It means that the fluctuations of a are higher for larger eddies 
with longer lifetimes - as it should. 

The experiment V works with largest correlation times and longest lifetimes, /corr = 8 and Tcorr = 
80. The birth rate i?tur is 0.0003 per time step - really low in comparison with the previous 



cases. Fig. g6| presents the time series of studied quantities. The resulting fluctuations of the alpha 
coefficients are extremely high. The mean value of r/x is 1.793 (4.5 • 10^^ cm^/s). The ratio 
SriT^ = 1.27 indicates the significant growth of the fluctuations compared with the mean value. 
The obtained ratios are 5q,^^ = 8.94 and Sayy = 6.62. The averaged in time values of axx and 
ayy decrease again to -0.015 (-1.5 • lO^cm^/s) and -0.017 (-1.7 • lO^cm^/s). 
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Figure 25: Time series for turbulence intensity, eddy diffusivity and alpha-tensor components 
and ayy for turbulence model A. For time, velocity (u and a) and diffusivity cgs units result 
after multiplication with 2.5 • 10^ s, 10^ cm/s and 2.5 • 10^"^ cm^/s. 
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Figure 26: The same as in Fig. p5| but for model V. 
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Table 3 summarizes the main results of our simulations. The fluctuations in the time series become 
more and more dominant with the decreasing number of eddies together with the increasing vortex 



size. The extreme case is plotted in Fig. |26|. For such models the fluctuations of both the eddy 
diffusivity as well as the a-effect are much higher than the averages. Even short-lived changes of 
the sign of the quantities were found. 

The a-effect fluctuations exceed those of the eddy diffusivity in all our models. The latter proves to 
be more stable than the a-effect against dilution of the turbulence. This is a confirmation, indeed, 
for those studies in which an a-effect time series is exclusively used in dynamo computations. 

8.3 A Plane-layer Dynamo 

Quite a few publications deal with the influence of stochastic a-fluctuations upon dynamo-gener- 
ated magnetic fields for various models. Choudhuri (1992) discussed a simple plane-wave dynamo 
basically in the linear regime. The fluctuations adopted are as weak as about 10 %. While in the 
ail -regime the oscillations are hardly influenced, the opposite is true for the a^ -dynamo. In the 
latter regime the solution suffers dramatic and chaotic changes even for rather weak disturbances. 
In the region between the both regimes the remaining irregular variations are suppressed in a model 
with nonlinear feedback formulated as a traditional a-quenching. 

A plane ID a^il-model similar to that of Section 6.1 is used to illustrate the influence of the 
fluctuating turbulence coefficients on a mean-field dynamo. It is not an overshoot dynamo (cf. 
Ossendrijver and Hoyng, 1996). The plane has infinite extent in the radial {x) and the azimuthal 
(y) direction, the boundaries are in the z-direction. Hence, we assume that the magnetic field 
components in both azimuthal direction, {B), and radial direction, {dA/dz), depend on z only. 
Note that z points opposite to the colatitude 6. 

The a-tensor has only one component (given in Section 2.3) vanishing at the equator, and it is 
assumed here to vanish also at the poles. The magnetic feedback is considered to be conventional 
a-quenching: a = aao^'(i?tot) with ^(Stot) oc B^^. The diffusivity is spatially uniform. The 
normalized dynamo equations are given in ( |62| ) and (|63|). 

Positive Ca describes a positive a in the northern hemisphere and a negative one in the southern 
hemisphere. (7q defines the amplitude of the differential rotation, positive Cq, represents positive 
shear dQ/dr. The dynamo operates with a positive a-effect in the northern hemisphere. Our 
dynamo numbers are = 5 and Cq, = 200. The turbulence models A, B and V are applied, 
flow patterns with small {A), medium {B) and with very large {V) eddies are used. 
A magnetic dipole field oscillating with a (normalized) period of 1.46 (Fig. ^7|a) is induced by a 
dynamo model without EMF-fluctuations {N — > oo, called O). The period corresponds to an ac- 
tivity cycle of 8 years in physical units. Also the turbulence model B produces an oscillating dipole 



but with a much more complicated temporal behaviour (Fig. 27c). It is not a single oscillation, the 



power spectrum forms a rather broad line with substructures. The 'quality'. 



Q = -P^, (74) 



AWcyc 



of this line (with Awcyc as its half-width) proves to be about 2.9. This value close to the observed 
quality of the solar cycle is produced here by a turbulence model with about 100 eddies along the 
equator (see also Ossendrijver et al, 1996). There are also remarkable variations of the magnetic 
cycle amplitude. 
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Figure 27: Dynamo-induced magnetic toroidal fields for the turbulence models O (27a), A (27b), 
B (27c) and V (27d). LEFT: time series, RIGHT: power spectra. The real time unit is 2.7 years. 
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The turbulence field V lead, not suprisingly, to a highly irregular temporal behaviour in the mag- 
netic quantities. Its power spectrum peaks at several periods (Fig. ^7|d). The overall shape of the 
spectrum, however, does not longer suggest any oscillations. The power of the lower frequencies 

—5/3 

is strongly increased, the high-frequency power decreases as Ucyc like a Kolmogoroff spectrum, 
indicating the existence of chaos. 

Our experiments basically show how the dilution of turbulence is able to transform a single-mode 
oscillation (for very many cells) to an oscillation with low quality (for moderate eddy population) 
and finally to a temporal behaviour close to chaos (for very few cells). It is possible that the 
observed low quality (33) of the solar cycle indicates the finite number of the giant cells driving 
the large-scale solar dynamo. Implications for the temporal evolution of the solar rotation law 
should be another output of such a cell number statistics. The consequences for rotation will be an 
independent test of the presented theory leading to the same number of eddies contributing to the 
turbulent angular momentum transport. For only 'occasional' turbulence a nontrivial time series 
for both the turbulence EMF-coefficients and the magnetic field is an unavoidable consequence. 
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